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ABSTRACT 

The shape of the primordial matter power spectrum encodes critical information on cosmological parameters. 
At large scales, in the linear regime, the observable galaxy power spectrum Pobs(^) is expected to follow the 
shape of the linear matter power spectrum P\ia{k), but on smaller scales the effects of nonlinearity and galaxy 
bias make the ratio Pobs(^)/^'iin(^) scale-dependent. We develop a method that can extend the dynamic range 
of the primordial matter power spectrum recovery, taking full advantage of precision measurements on quasi- 
linear scales, by incorporating additional constraints on the galaxy halo occupation distribution (HOD) from 
the projected galaxy correlation function Wp(rp). We devise an analytic model to calculate observable galaxy 
power spectrum Pobs(^) in real-space and redshift-space, given P\in{k) and HOD parameters, and we demonstrate 
its accuracy at the few percent level with tests against a suite of populated A^-body simulations. Once HOD 
parameters are determined by fitting Wp(rp) measurements for a given cosmological model, galaxy bias is 
completely specified, and our analytic model predicts both the shape and normalization of Pobs(fc). Applying 
our method to the main galaxy redshift samples from the Sloan Digital Sky Survey (SDSS), we find that the 
real-space galaxy power spectrum follows the shape of the nonlinear matter power spectrum at the 1 -2% 
level up to A: = 0.2/iMpc~' and that current observational uncertainties in HOD parameters leave only few 
percent uncertainties in our scale-dependent bias predictions up to A: = 0.5/iMpc~'. These uncertainties can be 
marginalized over in deriving cosmological parameter constraints, and they can be reduced by higher precision 
Wp(rp) measurements. When we apply our method to the SDSS luminous red galaxy (LRG) samples, we 
find that the linear bias approximation is accurate to 5% at A: < O.OS/zMpc"', but the strong scale-dependence 
of LRG bias prevents the use of linear theory at A: > 0.08/iMpc~'. Our HOD model prediction is in good 
agreement with the recent SDSS LRG power spectrum measurements at all measured scales {k < 0.2/!Mpc~'), 
naturally explaining the observed shape of Pobs(^) in the quasi-linear regime. The phenomenological "Q-model" 
prescription is a poor description of galaxy bias for the LRG samples, and it can lead to biased cosmological 
parameter estimates when measurements at A: > O.l/iMpc"' are included in the analysis. We quantify the 
potential bias and constraints on cosmological parameters that arise from applying linear theory and Q-model 
fitting, and we demonstrate the utility of HOD modeling of high precision measurements of Pobs(k) on quasi- 
linear scales, which will be obtainable from the final SDSS data set. 

Subject headings: cosmology: theory — dark matter — galaxies: halos — large-scale structure of universe 



1. INTRODUCTION 

In the linear regime, the power spectrum of matter fluctua- 
tions encodes information about the physics of early universe 
(e.g., the potential of the field that drives inflation) and about 
the matter and energy contents of the cosmos. The power 
spectrum of galaxies can be biased relative to the power spec- 
trum of matter (Kaiser 1984; Bardeen et al. 1986), but fairly 
general theoretical arguments imply that the shape of galaxy 
power spectrum should approach the shape of the linear mat- 
ter power spectrum Punik) at sufficiently large scales, i.e., 

PR{k) = blPUk)+No, (1) 



Harvard-Smithsonian Center for Astropliysics, Harvai'd University, 
60 Garden Street, Cambridge, MA 02138; jyoo@cfa.harvard.edu 
- Menzel Fellow 

^ Department of Astronomy, The Ohio State University, 140 West, 18th 
Avenue, Columbus, OH 43210; dhw@astronomy.ohio-state.edu 

^ Kavli Institute for Cosmological Physics, University of Chicago, 
Chicago, IL 60637; tinker@cfcp.uchicago.edu 

' School of Natural Sciences, Institute for Advanced Study, Einstein Drive, 
Princeton, NJ 08540; zhengz@ias.edu 

* John Bahcall Fellow 

' Theoretical Astrophysics Division, Los Alamos National Laboratory, 
Los Alamos, NM 87543; msw@lanl.gov 



where bo is a constant galaxy bias factor and 
PR(k) denote s the real-space galaxy power spectrurn 
("Coies "m?; 'Fry & Gaztanaga' 1994; 'Weinberg 1991 
Mann, Peacock & Heavens 1998; Scherrer & WeinberS 
1998; Narayanan et al. 2000; IBerlind & Weinberg! 120021 
iSchulz & White ,2006.) . The additive "shot noise" term 
A'o reflects both galaxy discr e teness and small scale 
clustering (Scherrer & Weinberg ' 1998t iMcDonaldl 120061: 
Smith, Scoccimarro & Sheth 2007); in general, it can differ 
from a simple Poisson sampling correction. In the linear 
regime, distortions of redshift-space structure by peculiar 
velocities also alter th e amplitude but not the shape of 
galaxy power spectrum (' Kaiseii[l987h . There have therefore 
been great efforts to measure the galaxy power spectrum 
on large scales from angular catalogs and redshift surveys 
and to use t he resu lts to tes t cosmo logical models (e.g., 
Yu & Peebl ej Il969t Baumgar t & Frvl [l 991; Feldma n et al. 
1994 iPark et al.1 fT994i: ,Linetal.l Il996t iSutherland et al. 
19991) . The enormous size of the Two Degree Field Galaxy 
Redshift Survey (2dFGRS; 'Colless et al.''2001) and the Sloan 
Digital Sky Survey (SDSS; York et al.. ,2000.) relative to 
earlier samples allows much higher precision measurements. 
The current state-of-the-art power spectrum measurements 
are Cole et a l.'s (2005) analysis of the power spectrum in 
the 2dFGRS, iPadmanabhan et al.l (l2007h . and Blake et al.'s 
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(2007) analyses of luminous red g alaxies (LRGs ) with 
photometric redshifts in the SDSS, and lPercival et alJ (l2007h 
and Tegmark et al.'s (2006) measurements from the SDSS 
redshift survey of main sample galaxies and LRGs. 

This paper investigates the problem of going from the 
galaxy power spectrum to the linear matter power spec- 
trum, and hence to cosmological conclusions. The latest 
observational analyses yield impressive statistical precision 
on scales near the transition from the linear to the non- 
linear regime, e.g., typical l-cr errors of 5-10% in P{k) 
at A; ~ O.lS/iMpc"'. The critical uncertainty in cosmolog- 
ical interpretation is therefore the accuracy of equation Q 
on these scales. The effects of nonlinearity and redshift- 
space distortions on the matter power spectrum can be com- 
put ed using num erical simulations or tuned analytic mod- 
els dSmith et alj |2003, and references therein), but details of 
galaxy formation physics can influence the rel ation between 
galaxy and matter power spectra in this regime. iPercival et alJ 
JIbO?) find that linear theory fits imply different cosmo- 
logical parameters if applied to measurements with k < 
0.06/iMpc~' or with k < O.lS/zMpc"', indicating that nonlin- 
ear e ffects have become significant in this regime. Further- 
more, 'Colilllnchir&^Sini (l2006h analyze the SDSS and 
2dFGRS galaxy samples and find that the measured shapes 
of galaxy power spectra differ at a level that cannot be ex- 
plained by the expected cosmic variance. They show that the 
likely source of the discrepancy is different scale-dependence 
of galaxy bias, originating from the different color distribu- 
tions of galaxies in the SDS S and 2dFG RS samples. 

Coleetal. ( 200 51), pTegmark et all (|2006|) and 
IPadmanabhan et al. l2007h approach this problem by 
fitting a parametrized model of scale-dependent bias, 

Ps-dk) = blPUk)^^^, (2) 

where we use Pgai(fe) to represent the galaxy power spectrum, 
which can be either in real-space or redshift-space. The func- 
tional form is devised for convenience to approximate the 
scale-dependent bias of galaxy sample s obtained by popu- 
lating the Hubble volume simulation (Evrard et al. 2002) us- 
ing a semi-analytic model of galaxy formation (Benson et al. 
120001) . Here A = 1.4/i-iMpc or l.T/j-'Mpc for real-space 
power spectrum P^ik) or angle-averaged redshift-space power 
spectrum Po{k) measurements, respectively, and Q is treated 
as a free parameter that is marginalized over in deriving cos- 
mological parameter constraints. This approach is adequate 
if equation (|2|i is a sufficiently accurate description of scale- 
dependent bias for some value of Q, but it could yield bi- 
ased parameter estimates or incorrect error bars if the actual 
scale-dependence is different. It also gives up on extracting 
cosmological information from scales where b ias mi ght be 
mildly scale-dependent. For example, iTegmark et al.l (l2006i 
hereafter, T06) find that cosmological parameters remain un- 
affected by changes in power spectrum measurements at A: > 
0.1/;Mpc~' once they marginalize over the value of Q. This 
implies that the statistical constraining power on cosmologi- 
cal parameters is lost at > 0. l/;Mpc~' by the marginalization 
process. 

In this paper, we present an alternative approach to re- 
covering the shape of the linear matter power spectrum, 
both more aggressive and more robust than "marginaliz- 
ing over 2" Our approach is based on the halo occu- 



pation distribution (HOD) framework, which describes the 
nonlinear relation between galaxies and matter by specify- 
ing the probability P(N\M) that a halo of mass M hosts 
number of galaxies of a given type, together with speci- 
fication of the relative spatial and velocity distributions of 
galaxies within halos.^ The HOD formalism has emerged 
as a powerful method of modeUng galaxy bias ( Jing et all 
^T99§, 'Seli akll2000[ |Ma & Frv"200o['Peacock & S mith.",200dr 
Scoccimarr o et al.l 2001; Berlind & Weinberg 2003) because 
the dynamics of dark matter halos can be accurately cal- 
culated using analytic approximations or A^-body simula- 
tions, and the effects of galaxy formation physics can be 
parametrized in terms of an HOD and inferred by fitting ob- 
servational data. 

Our strategy for extending recovery of the primordial mat- 
ter power spectrum is to use complementary information from 
the measurements of the projected correlation function Wp{rp) 
as a constraint to obtain HOD parameters given a cosmologi- 
cal model. We then predict the galaxy power spectrum Pg^iik) 
and study the scale-dependent bias 

b\k)=PgJk)/PM(k). (3) 

For each cosmological model, fitting Wp(rp) measurements 
determines HOD parameters and we can then compute a 
unique prediction of Psniik), both shape and normalization 
(which is essentially pinned to the amplitude of Wp(rp)). Un- 
certainties in HOD parameters introduce uncertainty in Pgai(fc) 
and b^{k), but these uncertainties can be accurately com- 
puted and marginalized over Therefore, we can extend the 
wavenumber range over which Pgai(fe) measurements can be 
used for cosmological parameter constraints, taking full ad- 
vantage of precision measurements on quasi-linear scales. In 
practice, we are just using the measured Pga\(k) and Wp(rp) to 
simultaneously constrain HOD parameters and the cosmolog- 
ical parameters and marginalizing over the former Relative to 
the 2-model approach, our method adopts a more physically 
motivated computation of Pgai(fc) and b^{k), requiring only the 
validity of the adopted HOD parametrization, and it brings in 
the additional information present in Wp{rp) rather than using 
only the P(k) shape itself to constrain the scale-dependence of 
bias. 

In principle, the power spectrum P{k) and correlation func- 
tion ^(r) contain the same information. However, they are 
in practice measured via different estimators and on different 
scales, where their signal-to-noise ratios are highest and sys- 
tematic errors are relatively well understood. The informa- 
tion in P(k) and ^(r) measurements on these non-overlapping 
scales is therefore not identical, but complementary. Fur- 
thermore, the projected correlation function Wp(rp) is mea- 
sured to ease the difficulty in interpreting nonlinear redshift- 
space distortion of correlation function measurements on 
small scales. Therefore, the addition of Wp{rp) measurements 
at rp < 30/!~'Mpc brings new information that is not present 
in P(k) measurements at A: < l/iMpc~^ 

A different approach to this problem is to develop an ana- 
lytic model for predicting the scale-dependence of galaxy bias 
by using higher-order perturbation theo ry (e.g., iMcDonaldl 
l2006t ISmith. Scoccimarro & Shethll2007i) . This approach is 

Throughout this paper, the term "halo" refers to a dark matter structure 
of overdensity p/pm — 200, in approximate dynamical equilibrium. 
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elegant and transparent in nature, since it is based on lin- 
ear theory and its extension to higher-order, while our ap- 
proach is less ab initio in the sense of incorporating elements 
calibrated by numerical A^-body simulations in our analytic 
model. However, the critical uncertainty for this approach 
based on higher-order perturbation theory is its applicability 
on quasi-hnear scales (> O.l/iMpc"'), where first-order linear 
theory is known to be inaccurate, but the measurement preci- 
sion is highest in practice. In contrast, our approach is fully 
nonlinear, and phenomenological in nature, so it can be ap- 
plied down to small scales, limited only by the point at which 
uncertainties in the HOD parameters introduce systematic un- 
certainty in the P{k) recovery. 

Analyses of galaxy redshift surveys typically estimate the 
angle-averaged power spectrum Poik), i.e., the monopole of 
the redshift-space p ower spectrum (e.g., ICole et alj 120051: 
iPercival et al.l 120071) . Redshift-space distortions do not al- 
ter the shape in linear theory, but they do chan ge the shape 
in the trans-linear regime (e.g.. lCole et al.|[r994l) . and finger- 
of-god (FoG) effects have impact out to large scales (e.g., 
Scoccimarro et al. 200 Ih . iPadmanabhan et all (200 7) and 
Blake et aT] ( 120071) deproject the angular clustering measure- 
ments of the SDSS LRG sample using photo-z catalogs to esti- 
mate the real-space power spectru n i, independ ent of redshift- 
space distortions. Tegm ark et al.l (l2004 1200 6) use a linear 
combination of the redshift-space monopole, quadrupole, and 
hexadecapole that recovers the real-space power spectrum in 
the linear regime. We will denote this "pseudo real-space" 
power spectrum Pz^R{k). The redshift-space power spectrum 
estimators can be applied directly to galaxy redshift data or 
applied after compressing FoG effects. We will investigate 
Pgai(fc) and b^{k) for all of these cases. 

To this end, we develop an analytic model in § |2] for cal- 
culating real-space and redshift-space galaxy power s pectra 
given P\iaik) and a galaxy HOD, drawing on the iTinke^ (l2007h 
model for redshift-space distortion, which i mproves on previ- 
ous work (e.g ., Seliak 2001; White 200ll lKang etani2002t 
ICooravl [20041) . i Tinker (2007) tests the model for comput- 
ing the redshift-space correlation function against a series of 
populated A^-body simulations. Here we extend the model 
and present additional tests of its applicability to modeling 
redshift-space power spectra in §[3] 

In this paper, we use HOD parameters for volume- 
limited galaxy samples that have well defined classes 
of galaxies, focusing on SDSS main galaxy samples 
with absol ute-magnitude limits Mr < -20 and Mr < -21 
(IZehavi et a l. 2005b) in § H and SDSS LRG samples with 
absolute-magn itude limits -23.2 < M„ ^ -21 .2 and -23.2 £ 
Mg < -21.8 (lEisenstein et al.1 1200 ll IZehavi et al.1 l2005at 
IZheng et"aLl 1200 8') in § |5] for application of our method.^ 
Mor e complete model ing of the conditional luminosity func- 
tion (lYang et al.ll2003h might allow use of flux-limited galaxy 
catalogs, though it requires more free parameters to provide 
complete descriptions of the galaxy samples. Here we only 
consider volume-limited galaxy samples, whose results can 
be combined to improve statistical precision. We summarize 
our main results in § |6] 

2. CALCULATIONAL METHODS 
2.1. Numerical Model 

' For brevity, we quote the absolute magnitude thresholds M, - 5 log h and 
Mg-5 log h for /i = 1 . 



We use the A^-body simulations of [Tinker et al.l (l2006h to 
test our analytic model calculations of the correlation func- 
tion and the power spectrum in real-space and redshift-space. 
These are five simulations of a flat A CDM universe usin g 
the publicly available tree-code gadget ( Springel et al.ll200Tl) . 
and all the simulations are performed with identical cos- 
mological parameters except for the random seed numbers 
used to generate initial conditions. The initial scale-invariant 
(» t =1) power spectrum is modified by the transfer function 
of Efst athiou et al.l d 19921) with shape parameter F = 0.2. The 
simulation was evolved from an expansion factor a = 0.01 to 
fl = 1.0 with n„ = 0.1, flA = 0.9 and erg = 0.95 at z = 0. To 
cover a range of parameter space spanned by ft,,, and crs, we 
use earlier outputs to represent different cosmological models 
from the simulations. Our choices for the earlier expansion 
factors are Oout =0.84, 0.64, 0.49, and 0.40. These outputs 
correspond respectively to simulations with different parame- 
ter combinations (a„, cr8)=(0.16, 0.90), (0.30, 0.80), (0.48, 
0.69), and (0.63, 0.60) with the identical power spectrum 
shape (F = 0.2) but evolved beginning at expansion factor 
a = 0.01 /flout- This procedure correctly provides the density 
field that would be obtained from an independent simulation 
evolved to z = w ith the correspond ing parameter combina- 
tion of fl,„ and (Tg (IZheng et al.ll2002|) . We evolve 360^ parti- 
cles in a volume of comoving side length 253 /i"'Mpc to take 
into consideration that the lowest mass halos that host galax- 
ies with Mr < -20 contain at least 32 particles. Dark matter 
halos are identified by using the friends-of-friends algorithm 
(FoF: lDavis et aljri985h with a linking length of 0.2 times the 
mean interparticle separation, i.e., 140/;"' kpc. 

To populate dark matter halos with galaxies in A^-body 
simulations, we use HOD parameters listed in Table [T] that 
are chosen to match the mean number density ilg and pro- 
jected correlation functions Wp{rp) of the SDSS galaxy sam- 
ples with absolute-magnitude limits Mr < -20 and Mr < -21 
JZehavi et al. 2005b). In our standard HOD parametrization, 
the number of central galaxies is a step function changing 
from zero to one at a minimum halo mass M^in. Therefore, 
halos of mass M < Mmin lack galaxies. We assume (A^sat) « M 
at high masses with a smooth cutoff at low mass. Therefore, 
the number of satellite galaxies is, 

for a halo of mass M > M^in, and the distribution of satel- 
lite galaxy nurn ber f (A^satl (A^sat)) is assumed to be Poisson 
jKravtsov et al.ll2004 IZheng et al.ll2005h . This parametriza- 
tion is well suited to our purposes, but we also investigate the 
effect of adopting a more flexible HOD parametrization in §|4l 
We replace ha los identified by the FoF algor ithm by spheri- 
cal NEW halos ( Navarro. Frenk & Whi"t3ll997i) with identical 
mass, truncated at virial radius Ryh, within which the mean 
density is 200 times the mean matter density. The concentra- 
tion parameters Cdm of dark matter halos are computed using 
the relation of iBuUock et al.l (120011) and are scaled to account 
for the different definition of halo overdensity adopted here. 
This NFW-replacement method reduces numerical artifacts 
caused by finite force resolution in our simulations. We place 
a central galaxy at the center of mass of halos. Assuming that 
satellite galaxies trace the dark matter distribution within ha- 
los, we place satellite galaxies following the NEW profile of 
halos. 
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In redshift-space, galaxies are displaced because of pecu- 
liar velocity. Central galaxies are assumed to be at rest rela- 
tive to the halo center; no velocity bias is assumed for central 
galaxies. For satellite galaxies, we add line-of-sight veloci- 
ties drawn from a Gaussian distribution with zero mean and 
dispersion 



/ GM 

V 2^vir 



1/2 



(5) 



to the velocity of the halo center of mass. This procedure is 
exact for isotropic singular isot hermal halos and is reasonably 
accurate for NFW profiles (see [Tinker et al.ll2006l for detailed 
tests). 

Finally, we compute the density contrast field by cloud-in- 
cell weighting the particle distribution onto 360'' grids and 
use the publicly available fast Fourier transform (FFT) code, 
FFTW, to obtain the Fourier components in units of the funda- 
mental mode of our simulation box, AA: = 0.02/iMpc~'. We 
deconvolve the cloud-in-cell weighting function Wcicik), and 
subtract shot-noise contributions 1 /Nan\ to compute the power 
spectrum at each k: 



(6) 



gal 



where the weighting function is 



Wcic(^) = 



n 

i=l 



sm(-Kki/2kM) 



-I 2 



(7) 



and A^gai is the number of galaxies, ki is the i-th component of 
wavenumber k, and k^ = 4.5/!Mpc~' is the Nyquist wavenum- 
ber of our simulations. For computations in redshift-space, 
we simply displace particles using the z-component of the pe- 
culiar velocity scaled by the Hubble constant as they would 
appear to a distant observer at z = -oo. This procedure sat- 
isfies the distant observer approximation we adopt here, and 
it ensures periodic radial velocity fields in the simulation vol- 
ume, appropriate for FFT. Redshift-space multipoles are ex- 
tracted by least-squares fitting to the Legendre polynomial co- 
efficients (eq. lTSl ). We repeat the procedure for the x- and y- 
axes, treating each axis as the line-of-sight, and we average 
the resulting power spectra over the three line-of-sight direc- 
tions. 

2.2. Analytic Model 

Our analytic calculation of t he real-space galax y auto- 
correlation function ^R{r) follows Tinke r et aTl (l2005l) . which 
improves the method of Zheng (2004) with more accurate 
treatments of scale-dependent halo bias and halo exclusion. 
For a given galaxy sample with its projected correlation func- 
tion measurements Wp{rp), we obtain HOD parameters by fit- 
ting the mean space density Hg and Wp{rp), computed by 



(8) 



Wp{rp) = 2 / t/z^R 



ir^p + z')"' 



(9) 



where dn/dM is the halo mass function of iJenkins et alj 
( 1200 lb and we use Zmax = 4 0/i~'Mpc as adopted in SP SS clus- 
tering measurements (e.g.. lZehavi et alj|2004l200 5a b). The 
real-space galaxy power spectrum Pr(^) is computed by tak- 
ing the Fourier transform of the correlation CR(r).'" 

We compute the redshift-space correlation function 
^(fW , r±) using the probability distribution of galaxy pairwise 
velocities f(v,,r), 



dvdl+ar)]f(v,,rl 



(10) 



the 



line-of- 

1 / 



where r± is the projected separation, r|| 
sight separation in redshift-space, and = 100 kms~'(r|| 
z)//!"'Mpc is the pairwise velocity of galaxies separated 
by r = (r^+z^)'/^. Equation ( fTOl i is called the streaming 
model and h as been used to measure the mean matter density 
jPeacock et al. 2001; Cole et al. 2005) and to model redshift- 
space correlations (White.2001;,Seljak.2()01i). It is valid in the 
linear and nonlinear regime "(iFisheiil 1 9951; IScoccimarrol2004 
provided that the correct f{v-,r) is used. We adopt the prob- 
ability distribution function of Tinker (2007) for the galaxy 
pairwise velocities, which is an analytic model with some el- 
ements calibrate d on A^-body s imulations. We refer the reader 
to the work by iTinkeii (l2007h for extensive discussion and 
tests. 

2.3. Redshift-Space Multipoles 

Using the analytic model, we compute the redshift-space 
correlation function i^(r||,r^) given P\\a{k) and a set of HOD 
parameters, and we expand ^(r||,rj^) with Legendre polyno- 
mials. 



(11) 



1=0 



where r = (r\ + '"y)'''^ and fi = ry /r is the direction cosine of 
the separation and line-of-sight vectors. The redshift-space 
multipole component is then 



21+1 



dfiLi{fj,)£,(r,fj,). 



(12) 



We use Li(jjl) to denote Legendre polynomials to avoid confu- 
sion with redshift-space multipole power spectra Pi(k) defined 
below. The redshift-space power spectrum P{k\^,k±) can be 
similarly expanded using Legendre polynomials. 

In general the correlation ^^(r) from galaxy pairs in two distinct halos 
is computed already in Fourier space, and one only needs to add the Fourier 
transform of the correlation £,^(r) from galaxy pairs in the same halo to com- 
pute PR(k). However, ^^(r) is modified to account for the finite extent of 
halos and the scale-dependence of halo bias, and hence the contributions to 
Pulk) from pairs in two distinct halos cannot be expressed in a si mple form 
in Fo urier space (for details, see iZheng||2004l : ITinker et alj|2005t lYoo et alj 
l200a . 
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p{k,fi)=Y,Liifi)m), (13) 

/=o 

where k = (k\ + k^^Y^^ and = k^^/k, in analogy to quantities 
in configuration space. The reflection symmetry of the cor- 
relation function and power spectrum ensures that multipoles 
with odd / vanish on average. Making use of the fact that 
^(^W , r±) and , k±) are Fourier counterparts, each redshift- 
space multipole component is computed by 

/•OO 

P,{k) = Am' r^^,ir)j,ikr)dr, (14) 
Jo 

where /"/(x) are spherical Bessel functions dCole et anil994t 
lHamiltonlll998h . Note that the quadrupole (/ = 2) components 
of Pi{k) and ^/(r) have opposite sign. 

Equation (fl4b requires knowledge of ^;(r) on large scales, 
while our analytic model is only tested at r < 40/!~'Mpc. 
Therefore, we compute ^;(r) at r > 40/i~'Mpc with the lin- 
ear approximation for redshift-space distortion. In the linear 
regime, the multipole expansion of (,ir^^,r±) has only three 
nonzero multipoles: monopole ^o, quadrupole ^2, and hexade- 
capole ^4 (.Kaiser . 1987 ). which are in turn related to ^Rir), 

6w=c2(eRW-ew), 

W) = C4 (CR('-) + 2.5e(r)-3.5|(r)) , (15) 

where Q = 1 + |/3+ C2 = 7/?^ C4 = (3 = 

^m^/^o, bo is the asymptotic galaxy bias factor, and the barred 
correlations are 

ar)=i / s^^R(s)ds, (16) 
Jo 



ar)=4 / *V*)«fs (17) 
Jo 

(lHamiltonlll992h . We compute the three multipoles at r > 
40/!~'Mpc, rather than ^(r||,rj^) itself. We first compute the 

values of Q at r = 40/;~'Mpc, where the deviations from the 
linear theory predictions are less than 5%, then we smoothly 
transition Ci{r) values to riin beyond which the adopted C/(r) 
values exactly become the linear theory predictions. The 
redshift-space multipoles ^;(r) at r > 40/!~'Mpc are then ob- 
tained by using equations ( fTSl l with C/(r), in place of con- 
stant C/, the linear theory prediction. We simply set r\i„ = 
500/j~'Mpc, and Pi{k) values are insensitive to t he choice of 
riin as long as r^^ > 100/;~'Mpc. We adopt the ISmith et al.l 
((2003") prescription for the matter power spectrum Pmmik) and 
its Fourier transform ^,„,„(r) to compute ^R(r) on large scales, 
and the asymptotic galaxy bias factor is computed by 

bo=^ f dM^(N)Mbk{M), (18) 
rig Jq aM 

where bh(M) is the bias factor of halos of mass M. We use the 
ISheth. Mo & Tormen (2001b) formulation with coefficients 
obtained by Tinker et al.i t2005i) . which yield a better fit to 
the simulations. 
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■ M_,S-20 sample 
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Tp (h^'Mpc) 

Fig. 1 . — Projected correlation functions Wp(rp) for SDSS galaxy samples 
with absolute-magnitude limits Mr < — 21 (upper pointsAines) and M, < —20 
(lower points/lines). Solid lines represent the best-fit analytic model predic- 
tions of Wpixp) for five different combinations of cosmological parameters, 
with increasing eg and decreasing Qm from lowest to highest curves (see Ta- 
ble[D. Shaded regions show the statistical uncertainty of the central cosmo- 
logical model only, computed from the error on the mean of five independent 
A'-body simulations. For comp arison, we plot the Wp(rp) measurements of 
the galaxy samples from Zehavi et al,| 42p05tf), with eiTor bars that are the 
diagonal elements of the covariance matrix. The model fits are acceptable 
when the full covariance matrix is considered. 



3. RECOVERING THE REAL-SPACE GALAXY POWER 
SPECTRUM 

Before turning to the bias between the galaxy power spec- 
trum and the linear matter power spectrum, we investigate 
how well the method used by Tegmark et al. (2004, 2006) re- 
covers the true real-space galaxy power spectrum. Along the 
way, we test the accuracy of our analytic model prediction for 
the redshift-space power spectrum against the results obtained 
from the A^-body galaxy catalogs described in § 12.11 Our basic 
approach to predicting the galaxy power spectra is that we first 
determine HOD parameters for an observed galaxy sample 
given a cosmological model, then calculate the galaxy power 
spectra using this inferred relation between galaxies and dark 
matter halos. We obtain HOD parameters by fitting the mean 
number densities Hg and projected correlation function mea- 
surements Wp{rp) of the SDSS My < -20 and Mr < -21 galaxy 
samples, taking into account the full covariance error matrix, 
estimat ed through jackknife resampling of the observational 
sample dZehavi et al.ll2005 bV 

Figure [U shows the projected correlation functions of the 
two SDSS galaxy samples, where the solid lines are the an- 
alytic model predictions for the five different cosmological 
models (listed in Table[T]l, obtained by fitting the lZehavi et alj 
(l2005bl) measurements shown as symbols. The error bars 
show only the diagonal elements of the covariance matrix of 
the Wp(rp) measurements, and the analytic model fits to the 
measurements are acceptable over a wide range of parame- 
ter combinations ($!,„, erg) when the full covariance matrix is 
considered, since the errors between data points are strongly 
correlated. However, there is strong degeneracy between the 
shape r and spectral index of the power spectrum, and the 
best-fit HOD parameters are insensitive to F and n,, as dis- 
cussed in § 14.21 below. Changing the adopted value of Zmax 
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TABLE 1 

HOD Parameters of the five A'-body models 



Model 




0"8 




M, < -20 






Mr < -21 




M™„(r'Mo) 


Mi(h-'MQ) 


Mc„,(/!-'Mq) 


M„i„(/!-iM0) 


Mi(/i-'M0) 


Mc„,(/r'MQ) 


1 


0.10 


0.95 


2.95 X 10" 


5.37 X 10'- 


1.40 X 10" 


1.64 X 1012 


2.05 X 10" 


2.89 X 10'3 


2 


0.16 


0.90 


4.80 X 10" 


8.03 X 10'- 


1.37 X 10'3 


2.60 X 10'^ 


3.06 X 10" 


3.47 X 10'3 


3 


0.30 


0.80 


9.33 X 10" 


1.32 X lO''' 


1.38 X 10'3 


4.83 X 10i2 


4.82 X 10" 


4.90 X 10'3 


4 


0.47 


0.69 


1.44 X 1012 


1.62 X 10'3 


1.78 X 1013 


7.15 X 1012 


5.69 X 1013 


6.53 X 10'3 


5 


0.63 


0.60 


1.89 X 10'2 


1.60 X W" 


2.46 X 10" 


8.91 X 1012 


5.65 X 1013 


6.85 X 10'3 



Note. — The HOD parameters of the five W-body modeLs are determined to reproduce the .same clustering w^trp) of the SDSS galaxy 
samples with M, < -20 and M,- < -21, and to match the number densities «^ = 5.74 X 10""^ (Z?"' Mpc)"" for the Mr < —20 sample and rig = 
1.17 X lO""* (Ir'Mpcy^ for the M, < -21 sample, respectively. 




k (hMpc-1) 

Fig. 2. — Dimensionless real-space and redshift-space multipole power 
spectra A^(k) for the HOD parameters appropriate to Mr < -20 galaxies. 
Various curves represent the analytic predictions of the central model for the 
corresponding galaxy power spectra indicated in the legend, and shaded re- 
gions show the statistical uncertainty on the real-space galaxy power spec- 
trum computed from the five A^- body simulations. The light, short dashed 
curve, labeled Pz^K(k) (see, eg. 1191 ), is a linear combination of redshift- 
space multipoles that reduces to P^(k) in linear regime; it is largely obscured 
by the solid curve. Thick dashed curves represent the linear (long) and non- 
linear (short) matter power spectra. Note that the quadrupole P2(k) crosses 
zero a.tk = 0.3/iMpc"' , and —Poik) is plotted at larger k. 



has little effect on the HOD parameters inferred by fitting 
Wp(rp), though it affects the values of these fits slightly. 
We use Zmux = 40/!~'Mpc for further analyses. The shaded re- 
gions represent the statistical uncertainty in the mean value 
of Wp{rp), computed from the dispersion among the five inde- 
pendent populated A^-body simulations, for the central model 
= 0.3, (Tg = 0.8; see Table [T]!. Our analytic model pre- 
dictions for Wp(rp) agree with the A^-body results within 5% 
fractional differences. 

Using the analytic model and assuming the central cosmo- 
logical model, we illustrate the dimensionless power spectra 
A^(k) = PP{k)/2T:^ of the < -20 galaxy sample in Fig- 
ure 12] where the thick solid line represents the real-space 
galaxy power spectrum Pr(A:) and the shaded region is the 
statistical uncertainty in the mean value of Pr(A:). The fi- 
nite box size of the simulations puts a limit on Fourier modes 
k > ^box = 27r/Lbox = 0.025/!Mpc~' that we can measure from 



the populated halo catalogs. We test our analytic model 
predictions for galaxy power spectra below. The two thick 
dashed lines are the matter power spectra of the assumed 
cosmological model; the long-dashed line is the linear mat- 
ter power spectrum Piin(k) and the short-dashed line is the 
nonlinear matter power spectrum Pa\{k), computed by using 
the ISmith et al.l (2003) prescription. On large scales, PR{k) 
has the some shape as P]m{k) with normalization differing by 
b^, the square of the large-scale galaxy bias factor, as pre- 
dicted by the linear bias approximation (eq.|[T|). However, 
on scales k > 0.2/!Mpc~', where Pni(^) departs from Pii„(k) or 
A|^(fc) ~ 1, PR(k) follows more closely the Pniik) shape than 
the P\in{k) shape. 

We also present the analytic model predictions for the 
redshift-space multipole power spectra Aj{k) = k^ Pi(k) J2tt^ . 
The dotted, long-dashed, and dot-dashed lines in Figure|2]rep- 
resent the redshift-space monopole Po{k), quadrupole Piik), 
and hexadecapole P^ik), respectively. On large scales, co- 
herent peculiar velocities produce redshift-space distortion, 
where the overdense regions shrink and the underdense re- 
gions inflate in redshift-space. Therefore, when averaged 
over angle, Po{k) is larger than Pr(A:). However, on small 
scales, nonlinear collapse and random motions in virialized 
objects stretch systems along the line-of-sight, giving rise to 
the Finger-of-God (FoG) effect, which inflates overdense re- 
gions and depresses their density contrast. Therefore, Po{k) is 
smaller than P^ik) on small scales. 

In the linear regime, the redshift-space power spectrum can 
be written as P{k,ij) = { \ + l3ii^)^PR(k), and the real-space 
power spectrum can be reconstructed by using a linear com- 
bination of the three redshift-space multipoles to remove the 
unknown variable /3, 

Pz^rW = PQ{k)-U^2{k)+^-P,{k), (19) 

which exactly reduces to P^ jk) if the linear theory approx- 
imation holds (lKaiser|[l987h . The upper panel of Figure [3] 
tests the ability of equation ( |T9] l to recover the true P^ik). 
Pz^Rik) recovers PR(k) to 5% at A: < 0.4/zMpc"', while it 
substantially underestimates PR(k) at k > 0.4/iMpc~'. This 
pseudo real-space power spectrum Pz^r(^) can be used to 
correct for the effect of redshift-space distortions and to es- 
timate Pr(^). a ^-by-^ application of equation iT% is, es- 
sentially, th e procedure ref e rred to as the disentanglement 
approach in iTegmark et al.l (12004 '2006), and our Pz^r(^) 
corresponds to Pggik) in their notation. Their alternative, 
modeling approach, is to construct two more power spec- 
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Fig. 3. — The upper panel plots pseudo real-.space galaxy power spectra 
and redshift-space monopoles with and without Finger-of-God (FoG) com- 
pression relative to the analytic model prediction for Pr(/:), assuming HOD 
parameters that fit the SDSS galaxy sample with M, < —20. Thresholds for 
the FoG compression are 400 kms"' and 750 km.s"' , as indicated in the leg- 
end. The redshift-space monopoles are scaled by a constant factor predicted 
by linear theory to match Pr(<:) at large scales. The bottom panel shows the 
fractional difference between the analytic model calculations and simulation 
results for the corresponding galaxy power spectra. Shaded regions repre- 
sent fractional statistical uncertainty on Pr(<:) from the W-body simulations 
caused by the finite simulation volume. The other power spectra have larger 
statistical uncertainties because of finite volume effects on redshift-space dis- 
tortions. Note that Pz^^(k) with and without FoG compression follow each 
other too closely on large scales to be separated from the thin solid line. 

tra, namely Pgv{k) and Pvv(k), from redshift-space multipole 
measuremen ts, then solve fo r PR{k) with the linear theory ap- 
proximation (iTegmark et al.l l2004, 2006). The former method 
gives a more robust approximation to Pr(A:) because of the 
cancellation of deviations in the multipoles from the linear 
theory approximation on nonlinear scales, but it yields larger 
statistical errors because it effectively marginalizes over un- 
certainties in the other two power spectra at each k. In con- 
trast, the modeling approach provides decorrelated estimates 
of PR{k) with smaller error bars, but it is subject to system- 
atic error in the linear theory approximation when model- 



ing the other two power spectra. ICole et al] (l2005l) use th e 
Fourier-based met hod of [ Percival. Verde & PeacocM (l2004l) . 
which extends the iFeldman et al.i (119941) method to include 
luminosity-dependent galaxy bias, computing the monopole 
Po{k) in our notation. We also take Po(k) as one of our princi- 
pal estimates of Pr(A:). 

A significant part of the nonlinearity in redshift-space dis- 
tortions arises from the random motion of galaxies in viri- 
alized objects, which can be identified by finding groups or 
clusters. The FoG effects can be removed by moving all 
galaxies in the same virialized object to the same redshift. In 
practice, clusters of galaxies are identified by applying the 
friends-of-friends algorithm to galaxy positions in redshift- 
space. Two galaxies are assumed to belong to the same clus- 
ter if the density windowed through an ellipse (usually sev- 
eral times longer in the radial than transverse direction) is 
higher than an overdensity threshold, and the radial disper- 
sion of member galaxy positions is set equal to the transverse 
dispersion t o compress the FoG effe ct, if the former exceeds 
the latter (Tegmark et al.ll2004ll2006f) . Here we assume a per- 
fect process of compressing the FoG effects; all the halos 
with velocity dispersion txh greater than a threshold are com- 
pressed by setting = 100 kms~', and none of the halos with 
(Th smaller than the threshold are affected. The pseudo real- 
space power spectrum Pz^^^^ik) with FoG compression thresh- 
old (Th = 750 km s"' agrees with Pr(A:) to 5% at < 0.5/iMpc"', 
and these halos are easily identifiable in practice. With a more 
aggressive threshold dh = 400 kms"', P^oo^(^) recovers P^ik) 
to 2% at A; < l/iMpc"'. We also compare the Po(k) shape to 
Pr(^) after scaling the constant factor Co in § 12.31 The scaled 
Poik) is depressed by 5% at = 0. l/iMpc"' compared to PRik). 
FoG compression helps at ~ O.l/iMpc"', but the difference 
between PQ°°(k) and PR(k) reaches 5% at /fe ~ 0.2/zMpc"' and 
10%at;t~0.3/iMpc-'. 

Having shown the agreement between the analytic model 
predictions and the A^-body results in Wp(rp), we now test 
the accuracy of the analytic model predictions for real-space 
and redshift-space power spectra against the A^-body simula- 
tions. The bottom panel shows the fractional difference in 
galaxy power spectra between the analytic model and the A^- 
body results, where the shaded region shows only the statis- 
tical uncertainty in the mean value of PRik) from the simula- 
tions. Note that the uncertainties on the other power spec- 
tra are larger and are not shown. The analytic model cal- 
culation of PRik) is accurate to better than a few percent at 
k > O.OS/iMpc"', while it is difficult to assess the statistical 
significance at ^ < 0.08/iMpc~', where the simulations only 
have few Fourier modes due to the finite box size. Since lin- 
ear theory should become accurate on large scales, it would 
be surprising if the analytic model became less accurate on 
this regime. Our analytic model also provides accurate pre- 
dictions for Poik) and Pz^r(^), both with and without FoG 
compression, at > O.l/iMpc"'. 

4. RECOVERING THE LINEAR MATTER POWER 
SPECTRUM 

We now turn to our principal results, the scale-dependent 
bias relation b^ik) = Pobsik) / Pimik) between observable galaxy 
power spectra Pobsik) and the linear matter power spectrum 
P\inik). As potentially observable power spectra, we con- 
sider PRik) (inferred from the angular clustering power spec- 
trum), Pz^Rik), and Poik), with varying levels of FoG com- 
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Fig. 4. — Scale-dependent bias functions calculated for tlie M,- < —20 galaxy sample. Panel (a) shows the real space P^ik), and panels (b) to (d) show Pz^R(k) 
and Po{k) with varying levels of FoG compression. In each panel, circles show the ratio of the nonlinear matter power spectrum to the linear matter power 
spectrum, and dotted lines show the Q-model curves that best reproduce our predicted h^(k). 



pression for the latter two. Here we use HOD constraints for 
the Mr < -20 and < -21 samples of IZehavi et"aLl (l2005bh 
based on their Wp{rp) measurements from SDSS Data Re- 
lease 2. These could be further improved with Wpirp) mea- 
surements from subsequent SDSS d ata and with constr aints 
from the group multiplicity function (iBerlind et al.l2006l) . We 
investigate the uncertainties in the power spectrum recovery 
associated with our HOD modeling and with variations of 
the assumed cosmological model. We u se a cmbfast trans- 
fer function (ISeljak & Zaldarriagal ll996') to compute P\ia(k) 
for a given cosmology, in place of the Efstath iou et alJ (Il992h 
parametrization used in §|3]to test our analytic model against 
the A^-body simulations (which used these initial conditions). 
Unless explicitly stated otherwise, we use a cosmological 
model with n,„ = 0.3, = 0.7, «, = l.O,h = 0.7, nbh- = 0.02, 
and cTg = 0.9 in this section. 

4.1. Scale-Dependent Bias of the Mr < —20 and Mr < — 2 1 
Galaxy Samples 

Figure |4] plots the scale-dependent bias functions b^{k)/bfi 
of the galaxy sample with M, < -20, in which bo is the asymp- 
totic galaxy bias factor computed via equation (fTSl l. FigureSh 
shows the real-space PR{k) (solid line) and the nonlinear mat- 
ter power spectrum Pni(k) (circles), which follow each other 



remarkably closely. A suppression of Pai(k) at A; ~ 0.1/jMpc~ 
relative to P\in(k) results from the nonlinear damping of linear 
perturbati ons. The physical origin of the Pni(k ) shape is dis- 
cussed bv lSmith. Scoccimarro & ShethI (120071) . Note that we 
compute the nonlinear matter power spectrum Pni(^) using the 
ISmith et aP (l2003h prescription, modified to utilize a cmbfast 
transfer function. P\in(k) retains the baryonic acoustic oscilla- 
tions (BAO) imprinted by sound waves in the baryon-photon 
plasma before recombination, but non linear evolution washes 
out th e oscillations at higher k (see, lEisenstein. Seo. White! 
l2007h . The bias functions in Figure |4] exhibit small oscilla- 
tions at > 0. l/iMpc~' because the smoother nonlinear power 
spectra are divided by a reference P\in(k) that retains the BAO 
features at their original strength. 

When normalized to the large-scale amplitude, the real- 
space power spectrum of M,- < -20 galaxies falls below 
Piin(;t) by ~ 5% at ~ O.l/iMpc"', cHmbs above by 5% 
at A: ~ 0.22/!Mpc~', and rises rapidly thereafter However, 
for this galaxy sample, the assumption that PR(k) traces the 
nonlinear matter power spectrum remains quite accurate, 
to 1% at < 0.2/iMpc"' and 7% at = 0.3/zMpc"'. The 
dotted curve shows the Q-model prediction of equation (|2|i 
with Q = 10.6 (/!"'Mpc)^, which gives a least-squares fit to 
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Fig. 5. — Same as Fig.|4] but for the Mr < -21 sample. 



our predicted bias curve over the range O.Ol/zMpc ' < < 
0.3/iMpc"'. The largest difference is at A; ~ 0.04/zMpc"', 
where the Q-model predicts a 5% deviation from Pua{k) while 
we find 3%. The Q-model accurately traces the predicted bias 
shape beyond k = 0. l/iMpc"' . However, if one normalized the 
curves to match at A: ~ 0.04/iMpc~', which might well hap- 
pen in practice because of the large statistical uncertainties at 
low k, then deviations in the bias shape would be smaller at 
k < 0.1/iMpc"' and larger at A: > 0.1/zMpc"'. 

Figure ^ plots Pz^r(A) and Poik) with no FoG compres- 
sion. The bias shape for Pz^r(^) is qualitatively similar to 
the true real-space P^ik), but it follows the Piin(fc) shape more 
closely than the Pniik) shape at A: < 0.05/iMpc"'. Pz^R(k) 
is poorly described by the Q-model prescription, with differ- 
ences of 6% at A ~ O.OS/iMpc"'. The bias shape for Po(k) 
is completely different from those for PR(k) and Pz^r(^)- 
By ;t = 0.1/iMpc"', Po{k)/bl is below PiUk) by 10% and 
Pni(k) by 5%, and it remains below Prm(k) by 10-15% to 
k = 0.5/iMpc"'. The best-fit g-model has Q = 2.2 (h-^Mpcf, 
and it has several percent differences from our Poik) predic- 
tion at A; < O.l/iMpc"'. 

Figures]?]: and|4j/ illustrate the effects of FoG compression 
on Pz-,Rik) and Po{k). Changes to Pz^R(k) are minor, though 
the agreement with the true PRik) is significantly improved at 



k > 0.25/;Mpc Changes to the monopole Po(k) are much 
more substantial, especially at A > 0.2/zMpc~'; suppression 
of Po{k) by the virial motions of satellites is much stronger 
than for Pz^R{k). However, the shape of Po(k) remains quite 
far from that of the true ^r(A), even for the aggressive FoG 
compression with dh > 400 kms~' . 

Figure |5] plots the scale-dependent bias functions of the 
galaxy sample with M,- < -21, in the same format as FigurelD 
Compared to the < -20 galaxy sample, this galaxy sam- 
ple has a large-scale bias factor higher by 20%, and hence 
a power spectrum amplitude higher by nearly 50%. How- 
ever, the shape of b^{k)/b\ is nearly identical for the two 
galaxy samples over a large dynamic range k < 0.3/iMpc~', 
as one can see by using the nonlinear matter power spectrum 
(circles) as a reference. In each case, the best-fit ^-model 
reproduces our predictions for Pr(A:), Pz^Rik), and P^ik) at 
roughly the 5% level for k < 0.3/iMpc"' . 

4.2. Impact of HOD Uncertainties 

Our method for calculating Pobs(^) and then b^ik) is that we 
first determine HOD parameters by fitting Wp{rp) measure- 
ments given a cosmological model. With perfect knowledge 
of HOD and cosmological parameters, it should be possible 
in principle to calculate Pobsik) exactly and fit the observed 
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Fig. 6. — Impact of HOD parameter uncertainties on bias shape. We compare our best-fit HOD models (solid) to ten HOD models (gray) for each galaxy 
sample that have Ax^ < 1 relative to the best-fit models, randomly chosen from the Monte Carlo Markov Chain. Left-hand panels show mean occupation 
functions and right-hand panels show h^(k). 



galaxy power spectrum to constrain Punik) up to high k. How- 
ever, uncertainties in the parameters result in uncertainty in 
our Pobs(^) calculations. Therefore, the dependence of Pobs(^) 
on our adopted fiducial HOD model will increase uncertainty 
in cosmological parameter estimation relative to measurement 
errors alone. 

To investigate these uncertainties, we adopt the 5-parameter 
HOD model of Zheng et al. (2005) instead of the 3-parameter 
model described in § 12.11 so that we do not underestimate 
HOD uncertainties because of an overly restrictive form. We 
determine best-fit parameters using this parametrization and 
the CMBFAST power spectrum. Then we generate a Markov 
Chain Monte Carlo based on the Metropolis-Hastings algo- 
rithm by fitting the Wp(rp) measurements with covariance ma- 
trix. In this parametrization, the mean occupation function for 
central galaxies is 



mean occupation function for satellite galaxies is 



(A^cen)A 



1-1- erf 



logM-logMjr 



(20) 



where erfix) is the error function, which corresponds to a 
Gaussian scatter of width gi pgM in log^' at fixed halo mass 
dZheng. Coil & Zehavill2007l) . The sharp cut-off of § ITT] is 
the limiting case of aiogM 0. The distribution of A^cen about 
the mean is a nearest-integer or Bernoulli distribution. The 



(A^sat) 



M '■ 



M-Mq 



(21) 



forM > Mo, and halos of M < Mo are devoid of satelUte galax- 
ies, so the slope of the satellite occupation is now free instead 
of fixed to one as in § 12.11 The number of satellites is Poisson- 
distributed about (A^sat)M. 

In Figure |6] the left-hand panels show the mean occupa- 
tion functions of ten HOD models for each sample that have 
Ax^ < 1 relative to the best-fit models, randomly selected 
from the Markov Chain. The HODs are tightly constrained 
by the Wpivp) data, especially for the M,- < -21 sample. The 
right-hand panels show b^{k) of the ten HOD models. Here we 
present the HOD dependence of Pob&iK) in terms of b^{k) since 
P\\n{k) is identical for all the models, and we only show b^{k) 
for P^ik), noting that FoG effects are similar for the fixed cos- 
mological parameters. Changing HOD parameters in the l-a 
allowed range alters b^(k) by up to 4% at k = 0.5/iMpc~' for 
the Mr < -20 sample, and by < 2% for k < 0.2/iMpc"'. In- 
duced errors for the Mr < -21 sample are even smaller, by 
about a factor of two, because of the smaller measurement 
errors in Wp(rp). Hence its impact on cosmological parame- 
ter constraints is small at A: < 0.2/iMpc~', and we discuss the 
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Fig. 7. — b^(k) for cosmological model variations. For each model in the legend, we adjust HOD parameters to obtain the minimum x^"fit to Wp{rp) 
observations. We then compute Pfi(k) and ti^{k) = P[i{k)/Pi[^(k). The dashed curve in the left panel shows b-(k) for the model with erg = 1.0 scaled to match the 
large-scale bias factor of the erg = 0.8 model. 



105 ^ (a) p^Ji^l 




100 r 

1.4 - 
.2 1.2 ~~ 
g 1 - 

0.8 

0.01 



(b) P,(k) 



0.1 0.01 

wavenumber k (hMpc"') 




Fig. 8. — Effect of cosmological parameter variations on the linear matter 
(left panel) and observable galaxy {right panel) power spectra. The bottom 
panels show the variations relative to the fiducial model power spectrum (erg = 
0.9, H, = 1.0, ft = 0.7). 

impact as we include measurements at higher A: in § |5] Note 
that even these small uncertainties in Pobsik) associated with 
HOD parameter uncertainties can be reduced using Wp(rp) 
measurements from larger SDSS samples, and by bringing in 
additional constraints such as the group multiplicity function. 
These uncertainties are unlikely to make a major contribution 
to overall error budget in Punik) recovery. 

Our model fitting assumes that the galaxy HOD is inde- 
pendent of large-scale environment. However, recent studies 
show that the clustering of galaxy mass halos has a substan- 
tial d ependence on halo formation time as well as halo mass 
(e.g. JGao et aLl2005l:lHarker et al.ll2005h. ope n ing th e door to 
environment-dependent HODs. Croton et alj ( l2005h use the 
Millennium simulation of lSpringel e t al. (20051) to investigate 
the effect of environment-dependent HODs on galaxy cluster- 
ing in semi-analytic galaxy models and find that environmen- 
tal effects change galaxy bias factors for luminosity thresh- 
olded samples by a f ew percent, but with virtually no scale- 
dependence (see also lZhu et al.ll2006h . Our HOD parameters 
are insensitive to the large-scale amplitude of Wp(rp), being 
driven mainly by shape of Wp{rp) in a regime where the con- 
tributions of galaxy pairs from single halos and two distinct 
halos are comparable. Therefore, we s uspect that environ - 
ment dependence at the level predicted bv lCroton et alj ( l2005h 



would have few percent impact on the predicted asymptotic 
galaxy bias factor bl but probably much smaller effect on the 
scale-depend ence of b^{k)Jb},. Furthermore, studies of void 
probabilities dTinker et al.ll2007 ') and color-separated correla- 
tion functions in "scra mbled" group catalogs (Blantonet aQ 
|2006|) suggest that the ICroton et alj ( l2005l) models overpre- 
dict the environmental dependence of galaxy HODs. Nev- 
ertheless, the environment-dependent issue merits further in- 
vestigation in future work, as the statistical precision of the 
measurements themselves sets better 

4.3. Sensitivity of Galaxy Power Spectrum to P\m{k) 

In approach to real observations, we will need to simultane- 
ously fit cosmological parameters and HOD parameters using 
Pobs(fc) and Wp{rp) as constraints and to marginalize over HOD 
parameters. The ability to constrain cosmological parameters 
depends on the measurement errors in Pobsik) and on the sen- 
sitivity of Pohs(k) to Pun(k). Here we investigate the sensitivity 
of Pobsik) to Piin(fc) given the constraints of Wp(rp). We discuss 
the impact on cosmological parameter estimation in §|5] 

In Figures |7] and [8j we illustrate the changes in Pobs(^) and 
b^{k) for cosmological model variations using the < -21 
sample. We adjust HOD parameters to fit Wp{rp) with each 
cosmological model shown in the legend. Figure [TJj plots 
b^(k) for models with erg = 0.8 and erg = 1.0. Note that 
for specified cosmology, fitting Wp(rp) completely determines 
both the shape and amplitude of b^(k). The thick dashed curve 
shows b^(k) for the erg = 1.0 model scaled to have the same 
large-scale bias factor as the erg = 0.8 model. The few percent 
difference in scaled b^{k) at A: 0.1 -0.2/iMp c~' results from 
the shape variation in the ISmith et al] ( l2003h Pni(k) prescrip- 
tion for different erg values. 

FigureslTf) and|7J^ plot two more model sequences with vari- 
ations in rts and h, which change the shape of Punik) in similar 
but not identical ways. Roughly speaking, fitting Wp{rp) fixes 
the Pohs(k) amplitude at ^ ~ 0.5/zMpc"' for the Mr < -21 sam- 
ple. The large-scale bias factors are therefore slightly differ- 
ent for distinct n., or h values. At A: < O.l/zMpc"', the shapes 
of b^(k) are nearly identical and nearly scale-independent, 
among the and h model sequences. On smaller scales, the 
b^{k) curves separate as HOD parameters adjust to try to pro- 
duce similar Wp{rp) from models with different matter power 
spectra. 

Figure[8]plots Punik) and PRik) in the upper panels and their 




r {h"'Mpc) 

Fig. 9. — Test of analytic model predictions for real-space ^aCr) and 
redshift-space monopole 5o('')- Here we take the HOD parameters for the 
luminous red galaxy (LRG) sample from Zheng et al. ( 2008) as input and 
populate one large volume A'-body simulation at z = 0. Symbols represent 
measurements of ^n(r) and i;()(r) from the simulation and curves represent 
analytic model predictions. The attached bottom panel shows the fractional 
differences between the analytic model calculations and the simulation re- 
sults. Statistical uncertainties are computed by jackknife resampling of the 
eight octants of the simulation box and shown as eiTor bars for i;n(r) and 
shaded regions for ^o{r). 



ratios relative to the fiducial model power spectra in the bot- 
tom panels, for the same model sequences considered in Fig- 
ure|2l Thus, this figure shows the degree to which a change in 
the linear power spectrum produces a detectable change in the 
observable power spectrum once we allow the galaxy HOD to 
vary in a way that reproduces Wp(rp). A reduction in erg to 
0.8 can be almost exactly compensated by a change in the 
HOD — in linear theory, of course, such a change could be 
exactly compensated by the linear bias factor bo- Linear bias 
would not change the shape of the observable power spec- 
trum, but we see in Figure[8]that the PR{k) curves for different 
lis and h converge towards a common shape atk> O.l/iMpc"', 
though they resemble Punik) at larger scales. Thus, nonlin- 
ear dynamical evolution and the freedom to introduce scale- 
dependent bias via HOD variation reduces the discriminatory 
power of Pobs(^)- However, with the HOD constrained by 
Wpirj,), the amplitude of P^ik) is offset by 5% or more for 
0.1/zMpc"' < A: < 0.5/!"'Mpc in the = 0.9 and h = 0.6 mod- 
els. The T06 main galaxy Pobs(^) measurements have typi- 
cal uncertainty of about 10% per data point, so this level of 
offset could provide significant additional sensitivity to cos- 
mological parameters if the measurements are extended to 
k = 0.5/;Mpc~' . The caveat is that environmental dependence 
of the HOD (not allowed in our models here) might be able to 
change the large-scale amplitude of fobs(^) relative to Wp{rp), 
erasing the predicted offset. The systematic uncertainty due 
to this effect could be removed by adding and marginalizing 
over a multiplicative normalization factor (essentially a "bias 
offset") so that only the shape of Pobs(^) provides constraints. 
This technique, and the impact of plausible levels of environ- 
mental variation, merits further investigation in future work. 



AL. 




0.01 0.1 
k (hMpc-i) 

Fig. 10. — Test of analytic model predictions for Pr{/c) and Po(k), in the 
same format as in Fig.|9] 



5. BIAS AND POWER SPECTRUM OF LUMINOUS RED 
GALAXIES 

While we have focused on the SDSS main galaxy sam- 
ples so far, the most powerful measurements in galaxy power 
spectrum come from the luminous red galaxy (LRG) sam- 
ples (Eisenstein et al. 2001) because of the large effective 
volume that the LRG samples probe. We consider the bias 
shape and power spectrum of the LRG samples separately 
from those of the main galaxy samples for several reasons. 
First, LRGs are physically distinct from SDSS main galax- 
ies: they are mainly central galaxies, occupying massive halos 
of M > 1O'^^/i~'M0, and they appear to have much stronger 
scale-dependent bias (e.g., T06 find higher Q values). Second, 
constraints on HOD parameters and systematic uncertainties 
in their values have not been as extensively investigated; here 
we draw on parameter constraints from the recent study of 
IZheng et"aD (l2008h . Third, few simulations with the neces- 
sary large volume and dynamic range are available for test- 
ing our analytic model predictions for the LRG samples. The 
tests of our analytic model in this regime that we present here 
should be regarded as a first step. However, our modeling al- 
lows us to understand how the scale-dependence of LRG bias 
may differ from that of the SDSS main galaxies. More thor- 
ough investigation of HOD parameters and their uncertainties 
is necessary before applying our method to LRG Pohs(k) mea- 
surements to infer cosmologic al parameters. Her e we simply 
adopt HOD parameters from Zheng et alj (l2008l) and ignore 
the im pact of u ncertainties in the parameters for the moment. 

Zhe ng et al.l (|2008) obtain the HOD pai-ameters of the LRG 
samples with absolute-magnitude limit -23.2 < Mg < -21.2 
and -23.2 < Mg < -21.8 by matching the projected corre- 
lation function Wp(rp) and mean space density tig = 9.7 x 
10"^(/i"'Mpc)"^ and 2.4 x 10"^(/!"'Mpc)"^ and by account- 
ing for the error covariance matrix taken frorn IZehavi et alJ 
(I2005al) (see Appendix B of ' Zheng et all l2008l for the HOD 
parameters). The mean redshift of the LRG samples is z = 
0.3, and the two LRG samples are essentially luminosity- 
thresholded samples. The Mg < -21.2 sample is dominated 
by fainter galaxies, so it probes a smaller volume than the 
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Fig. 11. — Scale-dep endent bias curves for P^ik), Pz^R(k), and Po(k) of the luminous red galaxy (LRG) sample with —23.2 < Mg < —21.2, in the same format 
as Fig.|4] Adopting the'Zheng et al.'f2009) HOD parameters, we use the analytic model to predict the scale-dependence of LRG bias at z = 0.3. T he dotted curve 
in e ach panel plots the Q-ii^odel prescription with Q = 20 that approximately fits our calculations (this value is somewhat smaller than used in lTegmark et al.l 
12000 , see the text for details). Note that P„i(/:) shown as circles is less biased relative to Piin(^) at z = 0.3 than at z = 0. 



Mg < -21.8 sample. 

We use a new large A^-body simulation (IWarren et al.ll2006h 
to test our analytic model predictions for the LRG samples. 
The simulation is p erformed with the Hashed Oct-Tree code 
dWarren & Salmori| [T993), evolving 1024^ particles in a vol- 
ume of comoving 1086/!~'Mpc on a side from z = 34 to the 
present, using a ACDM cosmology (ft,,, = 0.3, ft^ = 0.7, 
Qb = 0.046, h = 0.7, n.5 = 1, and erg = 0.9). Dark matter ha- 
los are identified using the friends-of-friends algorithm with 
bfof = 0.2. Since we only have the simulation output at z = 0, 
we cannot directly compute quantities of interest at the mean 
redshift of LRGs (z = 0.3), but we can test whether our ana- 
lytic model is accurate. 

Figure |9] plots the real-space and redshift-space monopole 
correlation functions measured from the A^-body simulation 
at z = 0, populated by using the same HOD parameters of 
the LRG sample obtained by fitting Wp{rp) measurements at 
z = 0.3. The evolution of cosmic structure from z = 0.3toz = 
increases the abundance of high mass halos, so with these pa- 
rameters we set higher fig than observed. Note, however, that 
we use the simulation just for the purpose of testing our ana- 
lytic model. The error bars are computed by jackknife resam- 
pling of the eight octants of the cube. Our analytic predictions 
are in good agreement with the populated A^-body simulation. 
The marginally significant discrepancy at r = 4/i"'Mpc could 



indicate that our halo exclusion treatment is inaccurate at the 
10% level, but it is difficult to assess the statistical significance 
of the discrepancy with only one simulation. FigurefTOlshows 
the analytic predictions for PR(k) and Po{k), where the sta- 
tistical uncertainties are computed by using the eight octants 
as independent measurements. The analytic model provides 
good approximations to Pr(A:) and P{)(k) at k < 0.3/iMpc~'. 
The shot-noise power spectrum dominates the measurements 
of PR(k) and Po(k) at k > 0.3/iMpc~', and the discrepancy at 
k ~ 0.3/;Mpc~' may indicate that our shot-noise subtraction 
scheme is imperfect. Nevertheless, more simulations are nec- 
essary to better quantify the statistical significance of the devi- 
ation at A: ~ 0.3/!Mpc~' . The A^-body test at z = demonstrates 
that our analytic model can be used to compute the correlation 
functions and power spectra in the LRG regime with reason- 
able accuracy, and we suspect that these predictions would 
remain accurate at z = 0.3. 

Figures [TT] and [T2]plot the bias shapes of the LRG samples, 
computed by using our analytic model and accounting for the 
mean redshift z = 0.3 of the LRG samples. We now use the 
cosmological parameters = 0.24, n, = 0.954, h = 0.73, and 
0-8 = 0.75, consistent with the WMAP3 results, and LRG HOD 
parameters for the same cosmology. The bias shapes are in 
marked contrast to those for the Mr < -2 1 and M,- < -20 sam- 
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Fig. 12.— Bias shapes Pr(<:), i^_R(j(:), and Po(<:) of tlie luminous red galaxy (LRG) sample with -23.2 < A/j < -21 
dotted curves are computed with Q = 15. 
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pies, which would closely follow the nonlinear matter power 
spectrum {circles)}^ This is mainly because the fraction of 
satellite galaxies is small in the LRG samples, approximately 
5%, and because their host halos are massive (> 10'-'/i~'Mq) 
and hence highly clustered. The bright LRG sample is more 
biased by 15% compared to the faint LRG sample, but it 
shows less scale-dependence at A: < 0. l/;Mpc~' . The low frac- 
tion of satellite galaxies also suppresses the redshift-space 
multipoles arising from virial motions of satellite galaxies in 
halos, and the Pz^r(^) and Pdik) shapes have little difference 
compared to the Pr(A:) shape. In the bottom panels of Fig- 
ures [TT| and [121 relatively small changes in the bias shapes 
arise between two thresholds of FoG compression, because 
both cases basically suppress all the halos that host two or 
more LRGs. 

The dotted curves plot the Q-modeX prescriptions with 
2 = 20 and 15 (/i"'Mpc)^ for the faint and bright samples, 
which approximately follow our calculations. These values 
are smaller than the best-fit value Q = 30.3 (/z"'Mpc)^ found 
by T06. Because of the relative statistical weights of the P(k) 
multipoles, the quantity Pggik) that they measure is interme- 



' ' We use the ISmith et"an 120031) prescription to compute the nonlinear 
matter power spectrum. P„\{k) at z = 0.3 is nearly unbiased relative to P|in(^) 
at k < 0. 1/i Mpc ~', while at z = it is biased as is shown in Figs.|4]and|5] In 
Figs. llll and ll2l we plot shorter ranges along the .r-axis than in Figs.|4]and[5] 
because of the stronger LRG scale-dependence. 



diate between our Pz^R(k) and P()(k), but closer to the latter 
(see the discussion in § [3] above). The scale-dependence of 
LRG bias that we predict from HOD modeling is therefore 
weaker than that inferred by T06 by fitting the power spec- 
trum. The discrepancy may reflect the difference between the 
"defogging" procedure used by T06 and the perfect FoG com- 
pression assumed here, though since we find that FoG com- 
pression has little impact this explanation would imply that 
the T06 method overcorrects FoGs. (More direct evidence for 
such overcorrection is presented by B. Reid et al. [in prepa- 
ration]). It should also be emphasized that no values of Q can 
fit our calculations with 5% accuracy up to = 0.2/!Mpc~'. 
However, we do confirm a basic result of T06 model fit- 
ting: scale-dependence of LRG bias becomes substantial at 
k > 0.1/iMpc~', and it is much stronger than that of less lumi- 
nous, main sample galaxies. 

Proper application of our method to observations will re- 
quire several key steps beyond the scope of our current in- 
vestigation. First and foremost is careful matching of galaxy 
samples analyzed for Pohs(k) measurements and HOD mod- 
eling. T06 use "flux-limited" LRG and main galaxy cata- 
logs from SDSS DR4 to obtain Pobs(^) measurements and 
constraints on the scale-dependent bias, while we use Zheng 
et al.'s (2008) "volume-limited" galaxy samples to constrain 
HOD parameters and predict Pobs(^)- Second, as mentioned 
earlier, we need more large-volume simulations at z > for 
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Fig. 13. — Predicted LRG power spectrum and the SDSS measurements. 
Points are tlie power spectrum measurements of tlie flux-limited LRG sam- 
ple without defogging, taken from iTegmark et al. ( 2006) (see the text for 
the difference in the analyzed samples). Using our analytic model, we plot 
Pgg(k) = 0.8Po(<:) - 0.07P2 (<(: ) +O.OO6P4 {k ) whic h closely matches the power 
spectrum estimate used in Tegmark et alj j2006l) . after statistical weights on 
multipole measurements are considered. The shaded gray bands show our 
predicted power spectrum of the volume-limited LRG samples with the lower 
bound for the LRG sample with -23.2 < Mg < -21.2 and the upper bound 
for the LRG sample with -23.2 < Mg < —21.8. The solid line represents 
Plinik) of our fiducial cosmological mo del and the dashed line is the Q-model 
prediction with Q = 30.3 quoted in Tegmark et al.' (2006) that best fits their 
"defogged" measurements. The vertical dotted lines represent the "nominal" 
nonlinear scale quoted bv iTegmark et al., ^2006), beyond which the measure- 
ments provided little leverage on cosmological parameters. 

testing the accuracy of our analytic model in a regime ade- 
quate for LRG clustering. Third, we need to more fully in- 
vestigate HOD uncertainties associated with analytic model 
fitting of Wp(rp) and with possible environmental dependence 
of the LRG HOD. 

With these caveats in mind, we present preliminary re- 
sults on comparison of our analytic model predictions to the 
T06 measurements of the LRG power spectrum in Figure [T3] 
Since it is hard to correctly account for the difference in our 
FoG compression and their defogging processes, we instead 
use the Pobs(fc) measurements without defogging (see Fig. 22 
in T06 for the no-defog data). Points with error bars show the 
SDSS measurements. 

In Figure [T3] the shaded gray bands show P^^{k) ~ 
0.8Po('t)-0.07P2('t) + 0.006P4(/t) for the two LRG samples, 
constructed by using our analytic model predictions of Po{k), 
P2ik), and PA{k), with the lower bound computed from the 
faint sample and the upper bound from the bright sample (see 
the Appendix of T06 for the empirical formula for Pggik)). 
Considering the difference in the analyzed galaxy samples, 
the shaded region is likely to encompass the true Pobs(^), and 
it is indeed in good agreement with the measurements. Note 
that our analytic model prediction allows no freedom in the 
asymptotic bias bo, and hence in the normalization of Pggik), 
once HOD parameters are pinned down by fitting Wpifp)}^ 

Of course, it is no surprise that power spectrum and Wp(rp) measure- 
ments imply similar amplitudes of galaxy clustering, but this agreement is 
nonetheless reassuring. 



This is in contrast to the linear or Q-model fit, in which bo can 
be arbitrarily adjusted to fit the data. The solid and dashed 
lines show bfiP\inik) of the fiducial cosmological model and 
the 2-model prescription with Q = 30.3 (/!"'Mpc)^ quoted in 
T06, which best fits their "defogged" measurements.'-' At 
k < 0.09/;Mpc~', all predictions are consistent with the lin- 
ear model (solid) at roughly the 10% level, justifying the ap- 
phcability of the linear model with the same accuracy. How- 
ever, significant scale-dependence of LRG bias prevents the 
use of the linear model at A: > O.l/zMpc"'. The vertical dot- 
ted lines show the nominal nonlinear scale quoted by T06; 
they fit g-models up to = 0.2/iMpc~' but note that data be- 
yond k = O.l/zMpc"' provided little leverage on cosmological 
parameters and mainly constrained Q. 

We further investigate this point by considering variations 
of the fiducial cosmological parameters in Figure[T4] For each 
cosmological model variation, we first re-fit Wp{rp) to obtain 

best-fit HOD parameters and then compute Pggik). The solid 
lines are computed by averaging Pggik) of the two LRG sam- 
ples with weight by their number density. Figure [T4J3 shows 
a sequence of Pggik) with varying cg values. The predictions 
are virtually identical, showing that the combination of Pggik) 
and Wpirp) has no direct constraining power on erg given the 
flexibility of our 5-parameter HOD model. Since the predic- 
tions are obtained by applying our analytic model of Pggik) to 
Zheng et al.'s (2008) HOD fits to Wpirp), both of which de- 
pend on (jg in a complex way, this constancy of Pggik) is a 
reassuring consistency check. 

In Figures [T4b and [T4l:. we consider sequences with vary- 
ing «s and h, respectively. Once we have fit Wpirp) for a 
specific cosmology, there is no freedom to adjust the shape 
or amplitude of the model Pggik), except within the obser- 
vational uncertainties on the HOD parameters. We see from 
Figures [T4b and [T4b that normalizing to Wpirp) effectively 
forces different models to agree at ~ 0.08/iMpc~'. The 
largest separation among the models comes at large scales, 
k < 0.02/iMpc~'. However, the models also differentiate no- 
ticeably on smaller scales, especially when is varied. In 
linear theory fits, these scales would be discarded because the 
model predictions are unreliable. In Q-model fits, the infor- 
mation they contain would largely be lost by marginalizing 
over Q. The varying models also predict different Wpirp) 
for rp > 20/!~'Mpc, but these scales carry little weight in the 
HOD fitting. 

These results suggest that fitting Pobsik) to small scales with 
HOD parameter constraints from Wpirp) could yield signifi- 
cantly improved constraints on cosmological parameters. To 
investigate this point, and to compare linear theory, Q-model, 
and full HOD-model fitting, we generate synthetic Pobs(^) 
measurements equal to the predictions of our analytic model 
for the fiducial cosmological parameters and the correspond- 
ing Zheng et al. (2008) HOD parameters. We assign observa- 
tional errors that are a factor of two smaller than those found 
by T06 at the same value of k, roughly approximating the im- 
provement that might come from the future full SDSS data 
set.'"* At > 0.2/!Mpc~' (where T06 did not present mea- 

Since T06 provide best-fit 2- values for defogged data only, the dashed 
line should be considered as a reference, rather than a comparison. Note that 
their defogged data would roughly follow the dashed line. 

Error bars from the Baryonic Oscillation Spectroscopic Survey (BOSS) 
of SDSS-111 will be much smaller still. 
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Fig. 14. — Effect of cosmological parameter variations on the predicted LRG power spectrum, in the same format as Fig. 1131 For each variation of our adopted 
fiducial cosmological model shown in the legend, we compute power spectra of the two volume-limited LRG samples and solid lines are the average power 
spectra weighted by their number density. The sequence of Pgg(k) with varying erg in Panel (a) is nearly identical, while the predictions of Pgg(k) in Panels (b) 
and (c) decrease atk = O.OlftMpc"' with increasing n, or h. 



surements), we assign fractional uncertainties equal to those 
atA; = 0.2/!Mpc"'. 

Figure [15] shows the synthetic measurement and its errors 
with our analytic model prediction {dotted) in the left panels 
and the cosmological parameter constraints inferred by fitting 
these data using the three different models in the right pan- 
els. For simplicity, we fixed the combinations fl„,h^ = 0.1272 
and ilbh^ = 0.0222, which are best constrained from cosmic 
microwave background measurements, and we varied n, and 
h as two free parameters, assuming a flat ACDM universe. 
In Figure [T5b , the solid and dashed lines show the best-fit 
linear model and Q-model using the measurements only at 

< /t„ax = 0.09/!Mpc-'. We marginalize over bo for linear 
theory and (bo,Q) for Q-model fitting. Both models can de- 
scribe the data reasonably well (the dashed line is largely ob- 
scured by the solid line). The contours and shaded ellipses 
in Figure [TSt' show the 1-a confidence levels (Ax^ = 2.3) for 
all three models, and the symbols represent the best-fit pa- 
rameters. To compute the parameter constraints of the HOD 
model, we use the Fisher matrix formalism and marginalize 
over a cosmological parameter ag and two sets of five HOD 
parameters for each LRG sample. Note that analytic model 
fitting of Wp(rp) assumes a erg value, which is largely degen- 
erate in Pobs(^) measurements. 

The HOD model (shaded region) provides tighter con- 
straints than linear theory (dashed contour), because we 
have additional information from Wp(rp) measurements. The 
best-fit parameters (circle) of the linear theory fit are only 
marginally consistent at 1-cr level with the true values (as- 
terisk), because the linear bias approximation is no longer an 
accurate description of LRG bias even at ^ < 0.09/iMpc~', 
given the small uncertainties in Pobs(^) assumed here. The 
2-model prescription (solid contour) yields relatively less bi- 
ased best-fit parameters (triangle). However, marginalizing 
over Q compromises statistical constraining power, and this 
fit has the largest uncertainties in parameter estimates among 
the three models. 

The bottom panels compare the three models with A:inax = 
0.2 and 0.4/iMpc~', demonstrating the utility of higher k mea- 
surements. The best-fit Pobs(^) (dashed line) of linear theory 
in Figure [TSb shows large deviations, and the best-fit parame- 
ters (circles) in Figure [TSl:/ are significantly different from the 
true values (out of range Wj > 1 . 1 for A:n,ax = 0.4/zMpc~'). The 



linear theory fit prefers higher to compensate the deficit 
in Pohs(k) on small scales, while the overall Pobs(k) shape 
and baryon wiggles keep the h value relatively unchanged. 
Clearly, linear theory is an invalid description of Pobs(k) for 
LRGs when fc^ax > O.OQ/iMpc"'. 

A similar but less extreme trend of bias in best-fit param- 
eters is found when the Q-model prescription is applied to 
Pobs(k) with A-max > 0.2/iMpc~^ The best-fit parameters of the 
Q-model fit are (1.0, 0.74), deviating from the true value and 
from the best-fit values with Amax = 0.09/iMpc~' at the l-a 
level. Higher n, is again favored, but a higher Q can par- 
tially balance the trend of higher n,. With fc^ax = 0.4/iMpc~', 
the best-fit parameters are (1.09, 0.74) and biased at the 3- 
<T level. The g-model prescription also develops substantial 
bias in best-fit parameters as kj„^^ increase. In contrast, while 
our HOD prediction is accurate by construction in this exper- 
iment, it makes full use of the variations of Pobs(k) on small 
scales, providing tighter constraints on cosmological parame- 
ters as small scale measurements are included. 

Since we have generated the synthetic data from the HOD 
model, we have adopted the most optimistic scenario in which 
the HOD model with correct parameters can give a perfect de- 
scription of the data. Nevertheless, Figure[T5]conveys two key 
points. Fkst, (2-model fitting to the LRG Pobs(k) may yield 
unreliable parameter estimates for kmux = 0.2/zMpc~' (and lin- 
ear theory may be problematic even at A:n,ax = 0.09/!Mpc~'). 
Second, increasing /tmax from 0.09/iMpc"' to 0.2/zMpc"' with 
HOD model fitting can significantly reduce statistical un- 
certainties in cosmological parameters, and increasing to 
0.4/!Mpc~' produces a modest further gain. 

6. CONCLUSIONS 

We have developed an analytic model to predict observ- 
able galaxy power spectra Pobs(k) for specified cosmologi- 
cal and galaxy HOD parameters, and we have verified its 
accuracy using A^-body simulations. As potentially observ- 
able power spectra Pobs(k), we have considered the real-space 
Pr(A;), the redshift-space monopole Po(k), and the pseudo real- 
space Pz^R(k), with varying levels of Finger-of-God (FoG) 
compression for the latter two. Once HOD parameters are de- 
termined by fitting the number density ilg and projected corre- 
lation function Wp(rp) of the observed SDSS galaxy samples, 
given a specified cosmological model, our analytic model can 
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Fig. 15. — Constraints on cosmological parameters obtainable from a future "final SDSS" LRG sample, using linear theory, Q-model, and HOD model fitting. 
Left panels show the synthetic data with error bars and best-fit power spectra for the three models as smaller scale measurements are included in the data fitting. 
Right panels plot the l-cr confidence levels (Ax' = 2.3) and best-fit cosmological parameters {symbol) infen'ed by fitting the three models to the synthetic data. 



be used to predict Pohs(k) measurements. 

The large-scale normalization of our predictions is also 
fixed in the process of fitting Wp{rp), providing a unique pre- 
diction for each combination of cosmological and HOD pa- 
rameters. In practice, one can simultaneously fit cosmolog- 
ical and HOD parameters using Pobaik) and Wp(rp) as con- 
straints, then marginalize over the HOD in deriving cosmo- 
logical parameters. By implementing a complete physical 
model of nonlinear galaxy bias and drawing on the additional 
information in Wp{rp), our method allows one to take full ad- 
vantage of precision measurements of Pobs(k) on quasi-linear 
scales (k = 0.1 -0.4/!Mpc~'), where linear theory or the phe- 
nomenological Q-model may be insufficiently accurate. Our 
main findings are as follows: 

1. Our analytic mode l for calc ulating Wpijp) follows the 
method described in Tinker et alj ( [2006) , with the improved 
treatment of the scale-dependent halo bi as and e l lips oidal halo 
exclusion corrections. Drawing on the i TinkeJ (120071) model 
for redshift-space distortion, the analytic model is extended 
to incorporate calculating real-space and redshift-space power 
spectra. We have tested its predictions for Wpi/p) and Pobs(^) 
against populated A^-body simulations spanning cosmolog- 
ical parameter range Vl„, =0.1 -0.63 and <7% = 0.6-0.95, 
with HOD parameters matched to represent two SDSS galaxy 
samples with absolute-flux limits My < -20 and M, < -21 



dZehavi et al.ll2005"bl) . The analytic model reproduces the nu- 
merical results of Wp(rp) to 5% or better, and the predictions 
of Pobaik) are consistent with the numerical results to 2% 
at ;t = 0.1 - l/iMpc"' and to 10% at k = 0.025 -O.l/iMpc"', 
though the finite box size of the simulations makes it diffi- 
cult to assess the statistical significance of differences on large 
scales. 

2. For the M, < -20 galaxy sample, the pseudo real-space 
power spectrum Pz^Rik) recovers the true PRik) to 2% at 
k < 0.2/iMpc~', while the deviation between P^ik) and the 
scaled monopole Poik) is already 10% at k = O.l/iMpc"'. 
However, the deviation of Pz^R(k) from Pr(A;) becomes sub- 
stantial atk> 0.3/iMpc~' . This deviation can be partly reme- 
died by FoG compression, which suppresses nonlinear be- 
havior of the redshift-space multipoles caused by the ran- 
dom motions of satellite galaxies within halos. With FoG 
compression threshold cTh = 750 kms"', P^^^ik) can recover 
PRik) to 5% atk< 0.45/jMpc"', and at higher k for P^°\(k). 
FoG compression also reduces nonlinearity of the monopole 
power spectrum, but P^^ik) can only achieve 10% accuracy 
at A: < 0.3 /?Mpc~^ We c onclude that the pseudo real-space 
method of lTegmark et alJ ([1998) is an effective tool for recov- 
ering the nonlinear real-space galaxy power spectrum from 
redshift-space measurements, especially if it is combined with 
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accurate FoG compression. 

3. The nonlinear matter power spectrum describes the 
nonlinear real-space galaxy power spectra to 1% at < 
0.2/!Mpc"' for the M, < -20 and Mr < -21 galaxy sam- 
ples, up to an overall bias factor hfy The shape of the scale- 
dependent bias function b^{k)/bl forPz^R(^) is qualitatively 
similar to PK{k) ?A k < 0.3/iMpc"', but the shape for Paik) 
is completely different over the entire range we consider 
here. FoG compression makes little difference to b^{k)/b\ 
for Px^^{k), but a large difference for Po(^)- For these SDSS 
main galaxy samples, the Q-model prescription traces our cal- 
culation of P-gik) relatively well at A: > 0.1/;Mpc~', but its 
shape on large scales differs, so it might induce some overall 
bias in cosmological parameters when fitted to Pobs(^) mea- 
surements that have large uncertainties sA k < 0.05/;Mpc~'. 
Similar trends but with larger discrepancy are found in com- 
parison to our Pz^Rik) and P{){k) calculations. 

4. Uncertainties in computing Pob&{k) in our method arise 
from observational uncertainties in the HOD parameters and 
from uncertainty in the adopted parametrization itself. We 
have examined these uncertainties by adopting a flexible HOD 
parametrization with freedom to explore a wider range of 
plausibl e halo occupa t ion fun ctions. For the Mr < -20 sample 
with the lZehavi et alj (l2005bh uncertainties in Wp(rp), the un- 
certainty in the predicted Pohs(k) is 2% at A: = 0.2/iMpc~', be- 
comes progressively smaller at lower k, and climbs up to 4% 
at k = 0.5/iMpc~'. The uncertainty is a factor of two smaller 
for the Mr < -21 sample, roughly the ratio of the fractional 
Wp{rp) measurement errors of the two samples. We have not 
investigated the uncertainties assoc iated with possib le envi- 
ronmental variations of the HOD ( Croton et alj|2005b . Based 
on work to date, we expect that such variations might lead 
to few percent uncertainties in the overall normalization pre- 
dicted for Pobs(^) after fitting Wp{rp), but that the impact on 
scale-dependence of b^{k)/b\ would be smaller 

5. Moving to the LRG regime, we have tested our ana- 
lytic model predictions against the z = output of a large vol- 
ume, 1024^ -particle A^-body simulation (Warren et al. 2006), 
populated based on Zheng et al.'s (2008) HOD fits to Wp(rp) 
for tw o volume limited SDSS LRG samples dZehavi et alj 
l2005ah . The analytic model predicts fo('") and ^R(r) to 5% 
or better over the range 0.1/!~'Mpc <r< 30/!~'Mpc, and the 
predictions for Poik) and Pr(A;) have similar accuracy over the 
range O.Ol/iMpc"' <k< 0.3/iMpc"'. 

6. For the LRG samples, the linear (scale-independent) 
bias approximation remains accurate at the 5% level Xo k = 
0.08/iMpc"' for the -23.2 < M^, < -21.2 sample and tok = 
0.1/jMpc"' for the -23.2 <Mg< -21.8 sample. There is lit- 
tle variation among Pz^R{k), Po{k), and P^ik), because LRGs 
are mainly central galaxies in massive halos, so random mo- 
tions of satellite galaxies have little impact. Similarly, FoG 
compression has only a small impact on b^{k)/b\ for these 
samples. Both samples show strong scale-dependence of bias 
sA k > O.l/zMpc"', much more than for main sample galax- 
ies.'^ If we fit b^{k)/b\ from our HOD models with the best 
2-model over the range k = 0.01 -0.2/zMpc~', the largest de- 
viation is 7%. 

This result, obtained by fitting HODs and computing Pab^k) witli our 
analytic model, confirms the result of T06 infen'ed by fitting PohAk) with Q- 
models. However, the actual scale-dependence we find for LRGs is somewhat 
weaker than that inferred by T06. 



7. We have presented a preliminary comparison of our ana- 
lytic model predictions to the T06 measurements of the LRG 
^obs(^), with no FoG compression. The difference between 
our volume-limited samples and the T06 flux-limited sample 
precludes a full quantitative assessment, but the qualitative 
agreement is remarkably good over the full range of the mea- 
surements, k = 0.01 -0.2/iMpc"' (Fig. [T3]l. Fits with differ- 
ent cosmological parameters differ on large scales and, to a 
smaller degree, at A: > 0.2/!Mpc~', indicating that measure- 
ments to smaller scales would provide additional discrimina- 
tory power. 

8. Looking to the future, we have generated synthetic 
Pobs(^) data from our analytic model with error bars half those 
of T06, then fit them to successively higher A:n,ax with linear 
theory, the Q-model, and the HOD model. Cosmological pa- 
rameters from linear theory fits are badly biased for k^^^^ > 
O.l/zMpc"', while for k^^^ = 0.09/zMpc"' they are biased at 
less than l-cr. Parameters from the Q-model are minimally 
biased for fc^ax = 0.09/iMpc~', biased by 1.2-ct for fc^ax = 
0.2/zMpc~', and biased by many-cr for A-^ax = 0.4/iMpc~'. 
Since the synthetic data are generated from the HOD model, 
the HOD parameter estimates are unbiased, and the error 
bars in cosmological parameters shrink steadily as A;max is in- 
creased from 0.1/zMpc"' to 0.2/iMpc"' to 0.4/iMpc"'. 

Results 7 and 8 are especially encouraging. Using only 
the HOD model and the information in Wp(rp), our method 
predicts exactly the scale-dependent bias for LRGs that is re- 
quired to transform the linear power spectrum from WMAP3 
into the SDSS galaxy power spectrum measured by T06. This 
is in contrast to g-model fitting, where a phenomenological 
parameter (motivated by simulation results but with no clear 
physical interpretation) is introduced specifically to account 
for the difference between the linear theory P(k) and the ob- 
served power spectrum. 

Despite the clear evidence that scale-dependent bias affects 
the LRG power spectrum beyond k = O.l/iMpc"', result 8 
shows that one can gain substantial additional leverage on 
cosmological parameters with HOD modeling of power spec- 
trum measurements up to A: = 0.2-0.4/iMpc~\ and possibly 
beyond. Realizing this opportunity will require several in- 
vestigations beyond those presented here. First, we will need 
more large volume simulations to test and, if necessary, re- 
fine our analytic model to the level of accuracy demanded 
by the final SDSS data set. Second, we must explore more 
thoroughly the uncertainties associated with the HOD fitting, 
including alternative parametrizations, the impact of velocity 
bias on redshift-space predictions, and the possible impact of 
environmental variations of the HOD. Given the growth of 
current and future galaxy surveys in depth and redshift, these 
investigations will be needed to go beyond linear theory. Pre- 
cise measurement of the primordial matter power spectrum 
will play a crucial role in constraining cosmological parame- 
ters and testing dark energy theories. 
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